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7.1 A BRIEF HISTORY OF TIME, SIZE, AND SHAPE 


The rules by which anatomical size and shape are generated have intrigued scientists 
for centuries. In 1638, Galileo suggested a mathematical relationship between pro- 
portional changes in the shape of bones as animals increase in size, which he argued 
was a functional necessity for weight bearing (1914). The formalism of Galileo, 
whereby, physical forces and mathematical laws became integrated with studies 
of size and shape in biology, was most conspicuously encapsulated over a hundred 
years ago in the 1917 monumental tome by D’Arcy Thompson entitled, On Growth 
and Form (Thompson 1917). In a breathtakingly comprehensive manner, Thompson 
synthesized the observations of numerous predecessors and contemporaries, and 
through countless examples built a theoretical and experimental framework for 
describing changes in morphology that persists to this day (Stern and Emlen 1999; 
Arthur 2006). 

An essential component of Thompson’s treatise was his system of Cartesian coor- 
dinates he employed to map the geometrical transformation of organs and organisms. 
Many other biologists in the 1920s and 1930s were motivated to address similar 
questions on size and shape in both the scholarly and popular literature (Gayon 
2000). In 1926, John Haldane wrote a topical essay entitled, On Being the Right Size, 
in which he stated that, “The most obvious differences between different animals 
are differences of size, but for some reason, the zoologists have paid singularly little 
attention to them....For every type of animal, there is a most convenient size, and 
a large change in size inevitably carries with it a change of form” (p. 1) (Haldane 
1926). This was just one of many topics during Haldane’s career for which he showed 
remarkable prescience; by pointing out how little attention had been given to size and 
shape previously, he in effect anticipated a whole discipline. 


A Hux Ley AND ALLOMETRY 


Soon thereafter, Haldane’s close friend Julian Huxley in his Problems of Relative 
Growth (which he dedicated to Thompson) expanded the discussion of size and 
shape to include mathematical representations of morphological transformations 
that arise over time, specifically during ontogeny and phylogeny (Huxley 1932). 
Along with Georges Teissier, Huxley symbolized relative growth with an algebraic 
power formula and introduced the term allometry to explain how changes in shape 
relate to changes in size (Huxley and Teissier 1936; Gayon 2000). A major moti- 
vation of Huxley, as well as many others who followed was to gain insight into 
the developmental (e.g., genetic and cellular) mechanisms generating allometric 
changes in proportion during evolution (Hersh 1934; von Bonin 1937; Lumer 1940; 
Needham and Lerner 1940; Anderson and Busch 1941; Lumer et al. 1942; Clark 
and Medawar 1945; Rensch 1948; Huxley 1950; Reeve 1950; Kermack and Haldane 
1950; Bertalanffy and Pirozynski 1952; Gould 1966, 1971; Lande 1979; Alberch 
et al. 1979; Atchley 1981; Gould 1981; Coppinger and Coppinger 1982; Shea 1983; 
Atchley et al. 1984; Riska and Atchley 1985; Shea 1985; Coppinger et al. 1987; 
Deacon 1990; Godfrey and Sutherland 1995a; Coppinger and Schneider 1995; Stern 
and Emlen 1999; Smith et al. 2015). 
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Haldane and Huxley viewed size and shape predominantly through the prism 
of genetics, which during that period was surpassing embryology as the arbiter of 
acceptable explanations for mechanisms controlling the evolution of morphology. 
This grew from the seeding of Mendel alongside Darwin, and the paradigm being 
cultivated vis-a-vis genes, mechanisms of inheritance, and mutations that affect 
morphology from geneticists such as William Bateson, Richard Goldschmidt, and 
others (Bateson 1894; Bateson and Mendel 1902; Robb 1935; Goldschmidt 1938, 
1940). Ten years after Problems of Relative Growth, Huxley published Evolution, 
the Modern Synthesis (Huxley 1942) in which he somewhat unintentionally helped 
push embryology out of the field of evolution for almost thirty years. This does not 
mean that embryologists were not thinking about evolution at the time or thereafter, 
but genetics ruled the roost due primarily to the robust and highly visible efforts of 
some former embryologists like Thomas Hunt Morgan (Morgan et al. 1915; Morgan 
1919) and mathematical geneticists such as Ronald Fisher, Haldane, Sewall Wright, 
and Theodosius Dobzhansky (Fisher 1930; Wright 1931; Sinnott and Dunn 1932; 
Haldane 1932; Dobzhansky 1937). 

Huxley was well-versed in embryology and evolution, given that his close col- 
league, Gavin de Beer had written Embryology and Evolution in 1930, and Huxley 
coauthored Elements of Experimental Embryology with de Beer in 1934 (de Beer 
1930; Huxley and de Beer 1934). Even though Huxley’s Modern Synthesis has been 
viewed as a nail in the coffin for evolutionary embryology, Huxley left some open- 
ings for development to play a role. He stated, “The course of Darwinian evolution 
is thus seen as determined (in varying degrees in different forms) not only by the 
type of selection, not only by the frequency of mutation, not only by the past history 
of the species, but also by the nature of the developmental effects of genes and of 
the ontogenetic process in general” (p. 555) (Huxley 1942). Likely, this inclusion of 
developmental growth was influenced by his own earlier studies (Huxley 1924, 1932) 
and by those who continued to work contemporaneously on allometry (Hersh 1934; 
Gregory 1934; Lumer 1940; Anderson and Busch 1941; Lumer and Schultz 1941; von 
Bonin 1937; Needham and Lerner 1940; Lumer et al. 1942), as well as other mecha- 
nisms of evolutionary embryology, especially those advanced by Walter Garstang 
(1922) and de Beer (1930). 


B pe BEER AND HETEROCHRONY 


Gavin de Beer not only recognized the importance of allometry but also devised a 
series of definitions and schemas relating time to size and shape that is arguably one 
of the most important contributions in the history of the field of evolutionary develop- 
mental biology. First and foremost, de Beer was an evolutionary embryologist (Ridley 
1985; Hall 2000a; Brigandt 2006). His work emphasized the significance of changes 
in the timing of developmental events, or heterochrony, in transforming the morphol- 
ogy of a descendant relative to an ancestor. Heterochrony was initially conceived by 
Ernst Haeckel and has been applied in various scenarios to link development and 
evolution (Kollmann 1885; Russell 1916; Bolk 1926; Garstang 1928; de Beer 1930; 
Dechambre 1949; Gould 1977; Hall 1984; McKinney 1988b; Klingenberg 1998; Smith 
2003; Keyte and Smith 2014). de Beer classified eight modes of evolution through 
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which ancestral and descendant ontogenies can differ (de Beer 1930). He provided 
examples for each type of heterochrony but argued that neoteny, defined as the reten- 
tion of juvenile features in the adult form, was the one that truly allowed for large, 
rapid phenotypic change and morphological diversification. Other significant evolu- 
tionary concepts that de Beer advanced in the context of embryology include clandes- 
tine evolution, homology, and evolutionary plasticity. Notably, de Beer dropped the 
word Evolution that was in his first edition book title and instead adopted Embryos 
and Ancestors for the 1940 and subsequent editions (de Beer 1940, 1954, 1958), which 
can be seen as a reflection of how much the field of population genetics, and not 
embryology, laid claim to the study of evolution during that era. 


C THe Hoty Trinity OF TIME, SIZE, AND SHAPE 


Nonetheless, the theory that changes to the rate of growth and/or timing of events 
during ontogeny could alter the course of phylogeny continued as a subplot to the 
main story of evolution until becoming more generally accepted during the rebirth 
of evo-devo in the 1970s. Even Darwin in his Origin of Species was vexed and 
tantalized by the correlations of growth observed in embryos, which he acknowl- 
edged were a potential source of evolutionary variation (Darwin 1859). In Chapter [, 
Variation Under Domestication, Darwin wrote, “There are many laws regulating 
variation, some few of which can be dimly seen...I will here only allude to what 
may be called correlation of growth. Any change in the embryo or larva will almost 
certainly entail changes in the mature animal” (p. 11). He also stated: “If man goes 
on selecting, and thus augmenting, any peculiarity, he will almost certainly uncon- 
sciously modify other parts of the structure, owing to the mysterious laws of the 
correlation of growth. The result of the various, quite unknown, or dimly seen laws 
of variation is infinitely complex and diversified” (p. 12). Then again in Chapter V, 
Laws of Variation, Darwin explained that: “Changes of structure at an early age will 
generally affect parts subsequently developed; and there are very many other cor- 
relations of growth, the nature of which we are utterly unable to understand” (p. 168). 
Clearly, such correlations of growth were exactly on what Thompson focused, and 
his efforts helped lay the groundwork for a broad range of studies comparing changes 
in size and shape during development. 

All the more so, about a decade before Thompson’s seminal work, Charles Minot 
provided a complimentary and in many ways equally important embryological per- 
spective that connected the size of animals and/or their organs with the regulation of 
cell number, differentiation, and rates of growth as a function of age (Minot 1908). 
Borrowing from this idea, In Chapter III, The Rate of Growth, Thompson equated 
age with time and stated that “the form of an organism is determined by its rate 
of growth in various directions; hence rate of growth deserves to be studied as a 
necessary preliminary to the theoretical study of form, and organic form itself is 
found, mathematically speaking, to be a function of time” (p. 79) (Thompson 1952). 
Similarly, Huxley latched on to the importance of time when contemplating evolu- 
tionary changes in relative size. He proposed potential genetic mechanisms involving 
“(a) mutations affecting the primary gradient of the early embryo, on which the time- 
relations of antero-posterior differentiation depend; (b) mutations affecting specific 
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rate-genes; (c) mutations affecting specific ‘time-genes’—genes controlling time of 
onset and not rate of processes” (p. 242) (Huxley 1932). Accordingly, changes to 
these rate-genes and time-genes can affect growth gradients and alter morphology at 
multiple levels in a coordinated way. Such theories were supported by Goldschmidt’s 
discovery of genes that alter rates of development (Goldschmidt 1938, 1940; Dietrich 
2000), something which was also integrated into evolutionary embryology, and more 
specifically heterochrony, by de Beer (Hall 2000a). On this point, de Beer stated, 
“By acting at different rates, the genes can alter the time at which certain structures 
appear” (p. 20) (de Beer 1954). 

Therefore, primarily through the critical contributions of Minot, Thompson, 
Huxley, and de Beer during the first half of the twentieth century, the three param- 
eters of time, size, and shape became unified in essence as the holy trinity of evolu- 
tionary morphology. But while some of these authors and others strived to integrate 
findings from the emerging field of developmental genetics led by classically trained 
embryologists and morphologists such as Goldschmidt (1938, 1940, 1953); Conrad 
Waddington (1939, 1940, 1957b, 1962) and Ivan Schmalhausen (1949), a deeper 
understanding of the molecular and cellular mechanisms that unite time, size, and 
shape during ontogeny and phylogeny would have to wait for almost fifty years. 
Moreover, the neo-Darwinians remained very skeptical that developmental genetics 
could contribute to evolutionary theory, and thought-leaders such as Dobzhansky 
(1937, 1951) and Ernst Mayr (1963, 1983) argued most vociferously that all evolu- 
tion was microevolution arising from “the continuous adjustment of an integrated 
gene complex to a changing environment” (p. 332) (Mayr 1963). In other words, this 
was the prevailing synthetic theory that embraced natural selection and survival of 
the fittest, distribution of alleles at the level of populations, and gradual adaptive 
evolution as the sole agent of change. In this regard, the neo-Darwinians thoroughly 
rejected and even mocked the ideas of Goldschmidt (Gould 1982b), especially that 
small genetic changes affecting developmental time or rates could rapidly generate 
large phenotypic transformations in size and shape. 

So, by the 1950s, evolutionary studies predicated on allometry and heterochrony 
were either vastly overshadowed by the neo-Darwinian paradigm, or more pointedly 
they were viewed as gross oversimplifications of embryonic growth by developmen- 
tal biologists. Waddington, for example in a paper discussing how to measure size 
and shape in a meaningful and biologically relevant way stated that, “The validity 
of any biological conclusions which may be drawn from measurements of size and 
form depends far more on the adequacy of the physiological insight on which they 
are based than on the precision of the mathematical techniques used to summarize 
and compare them” (p. 515) (Waddington 1950). This sentiment begged the question 
of what governs growth over time and demanded a more in-depth probing of devel- 
opmental mechanisms regulating size and shape. 


7.2 TIME, SIZE, AND SHAPE REDUX 


Despite Waddington’s emphasis on acquiring a deeper understanding of developmen- 
tal processes and his admonishment of expending too much energy on generating more 
sophisticated and precise methods for measuring size and shape (Waddington 1950), 
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studies on allometry continued unabated for decades (Stern and Emlen 1999; Gayon 
2000). Moreover, a whole field of morphometrics burgeoned based on multivariate 
methods and ultimately computer-based algorithms for quantifying and visualizing 
complex changes in size and shape (Bookstein 1978, 1990; Benson et al. 1982; Siegel 
and Benson 1982; Marcus 1996; Zelditch 2004; Hallgrimsson et al. 2015). Granted, 
the technical ability to analyze size and shape became more refined over time, but 
results generally remained phenomenological. Therefore, many morphometricians 
endeavored to frame their studies within the context of quantitative genetics and/or 
heterochrony, in order to make predictions about mechanisms through which size 
and shape can change during ontogeny and phylogeny (Gould 1966, 1981; Lande 
1979; Alberch et al. 1979; Atchley 1981; Cheverud 1982; Benson et al. 1982; Riska 
1986; McKinney 1988a; Atchley and Hall 1991; Coppinger and Schneider 1995; 
Klingenberg 1998; Roth and Mercer 2000; Drake 2011; Smith et al. 2015; Lord 
et al. 2016). 


A CLOCKS FOR TIME, SIZE, AND SHAPE 


Some of the most prominent work, applying numerical methods to characterize growth- 
related changes in size and shape came at the end of the 1960s and in the 1970s from 
Stephen Jay Gould, who almost single-handedly made allometry and heterochrony 
fashionable again and also acceptable as alternatives to the adaptationist program for 
studying evolution offered by the neo-Darwinians (Gould 1966, 1971, 1977; Gould 
and Lewontin 1979; Gould 1981, 1982a; Gayon 2000; De Renzi 2009). 

Through a series of monographs and major papers on evolutionary allometry, 
Gould began to put developmental mechanisms front and center. Then, in a landmark 
book, Ontogeny, and Phylogeny, Gould (1977) traced the history of conceptual 
advances in understanding the relationship of development to evolution. Gould 
opened with the Great Chain of Being from the Greeks; continued to theories that 
ontogeny parallels or recapitulates phylogeny from various French and German tran- 
scendentalists such as Johann Meckel, Etienne Serres, Lorenz Oken, Louis Agassiz, 
and Ernst Haeckel; and finally described the outright rejection of recapitulation by 
Karl Ernst von Baer, Garstang, de Beer, and others. In the second half of his book, 
Gould revisited and expanded upon de Beer’s schema for heterochrony and pre- 
sented his own semi-quantitative clock model in which the hands for size and shape 
depicted the morphology of a species relative to its ancestor. Each clock allowed for 
size, shape, and age (i.e., time) to be altered separately during evolution and accord- 
ingly could be adjusted to represent the many manifestations of heterochrony such 
as neoteny, progenesis, pedomorphosis, proportional dwarfism, and proportional 
gigantism. 

In his impressive treatise and throughout his career, Gould confronted the neo- 
Darwinian view of morphological evolution head-on and argued forcefully for the 
role of development in macroevolutionary change (Gould 1966, 1971, 1977, 1982a, 
2002; Eldredge and Gould 1972; Gould and Lewontin 1979; Gould and Vrba 1982). 
But Gould was a paleontologist, not an embryologist, and one critical issue was 
that his clock model was essentially qualitative and static (like the models of de 
Beer), and defined simple evolutionary patterns or end states rather than capture the 
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complex and dynamic nature of developmental processes (Etxeberria and De la Rosa 
2009). Shortly thereafter, the embryologist David Wake and his 24-year-old graduate 
student Pere Alberch invited Gould to collaborate on what was to become an espe- 
cially celebrated paper that effectively launched the modern field of evolutionary 
developmental biology (Wake 1998; De Renzi 2009). 


B = QUANTITATIVE METHODS FOR TIME, SIZE, AND SHAPE 


Alberch thought Gould’s clock models were a good start in theory but insufficient in 
actuality (Reiss et al. 2008), and so in a paper entitled Size and shape in ontogeny 
and phylogeny, Alberch, Gould, Oster, and Wake (1979) presented a tangible quanti- 
tative method to describe the relationship between heterochrony and evolution. Their 
intention was to “clothe” Gould’s clock model in mathematics (Wake 1998) and in 
so doing build a better graphical framework for integrating changes in time with 
changes in size and shape during development and evolution. They formulated dif- 
ferential equations as a way to encapsulate a more dynamic view of heterochrony, 
which they described as shifts in the onset, cessation, or rate of growth, rather than 
as an end result (Etxeberria and De la Rosa 2009). 

This highly cited work became an instant classic that helped spawn a decade of 
conferences and books on how to measure size, shape, and time in the context of het- 
erochrony (Bonner 1982; Maderson et al. 1982; Raff and Kaufman 1983; McKinney 
1988b; Wake and Roth 1989; De Renzi et al. 1999; Reiss et al. 2008). Moreover, as 
part of the re-birth of evo-devo as a discipline, heterochrony became the lens through 
which all kinds of biology was viewed (Alberch 1980a; Alberch and Alberch 1981; 
Balon 1981; Coppinger and Coppinger 1982; Gould 1982a; Haluska and Alberch 
1983; Shea 1983; Bemis 1984; Hanken and Hall 1984; Roth 1984; Coppinger et al. 
1987; Geist 1987; Hoberg 1987; Slatkin 1987; Foster and Kaesler 1988; Hafner and 
Hafner 1988; Coppinger and Smith 1989; Roth and Wake 1989; Shea 1989; Coppinger 
and Feinstein 1991; Blanco and Alberch 1992; Zelditch et al. 1992; Blackstone and 
Buss 1993; Klingenberg and Spence 1993; Allmon 1994; Duboule 1994; Coppinger 
and Schneider 1995; Godfrey and Sutherland 1995a, b; Richardson 1995; Gilbert 
et al. 1996; Maunz and German 1996; Richardson et al. 1997; Smith 1997; Nunn 
and Smith 1998; MacDonald and Hall 2001; Vaglia and Smith 2003; Crumly and 
Sanchez-Villagra 2004; Tokita et al. 2007; Drake 2011; Nagai et al. 2011; Mitgutsch 
et al. 2011). 


C CONSTRUCTION RULES FOR TIME, SIZE, AND SHAPE 


While such morphometric approaches helped elucidate critical developmental stages 
and events whereby changes in size and shape occur, by necessity they often reduced 
the complex dynamic nature of development into something much more simplistic 
and static, their framework was typically applied globally at the level of organisms 
rather than in relation to individual systems or structures, they tended to divide con- 
tinuous development into artificially discrete steps in order to compare ontogenetic 
trajectories, and also, they largely left much to be understood in terms of underlying 
molecular and cellular mechanisms. 
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Seemingly anticipating these points and echoing Waddington’s sentiments, 
Alberch and his colleagues (1979) challenged the field when they expressed that: 
“We hope that our attempts to construct a quantitative theory will stimulate others 
to delve more deeply below the level of pure phenomenology and come to grips 
with the central issue underlying evolutionary diversification of size and shape— 
that is, the morphogenetic unfolding of genetic programs in ontogeny and their 
alteration in the course of phyletic evolution” (p. 297). 

Such an emphasis on the mechanistic and more dynamic aspects of development 
grew directly out of Waddington’s epigenetic landscapes and concepts like canaliza- 
tion, which basically served as metaphors for how gene regulation could alter the 
course of ontogeny and phylogeny (Waddington 1957a), and guided the remainder of 
Alberch’s remarkably influential but tragically foreshortened career (De Renzi et al. 
1999; Wake 1998; Reiss et al. 2008). To this very point, in his elegant first solo paper 
on the role of ontogeny in morphological diversification, Alberch (1980b) argued that 
“epigenetic interactions drastically constrain the universe of possible morphological 
novelties and impose directionality in morphological transformations through phy- 
logeny” (p. 654). In other words, even if a genetic mutation is random, the morpho- 
logical outcome is not. Why? Because developmental systems are highly integrated, 
iterative, accommodative, hierarchical, and ultimately defined by an “internal struc- 
ture” that limits “the realm of possible morphologies” (Alberch 1982a:319). 

In his subsequent and quite a formidable body of work, Alberch addressed the 
role of development in the evolution of size, shape, and other aspects of morphology 
on multiple levels in a wide range of organisms and organs. A critical concept that 
he advanced pertained to construction rules through which developmental systems 
are built and become altered from one morphological state to another during evolu- 
tion (Alberch 1982a, 1985; Oster and Alberch 1982; Oster et al. 1988). Accordingly, 
development consists of interwoven “dynamical systems, where a small set of simple 
rules of cellular and physicochemical interactions” lead to complex morphology 
(Oster and Alberch 1982:455). Evolutionary changes in organ size, for example, can 
be achieved by varying the quantitative parameters of cells, including the number 
of progenitors, the rate of proliferation, length of the cell cycle, and timing of dif- 
ferentiation. Other parameter values that can potentially be modulated pertain to 
“biochemical, cell-cell, or tissue interactions” (Alberch 1985:50), which, in turn, can 
affect developmental processes such as “rates of diffusion, mitotic rate, cell adhe- 
sion, etc.” (Alberch 1989:27). 

Throughout his research program, Alberch combined insights from comparative 
morphology, experimental embryology, and teratology, to generate models and other 
mathematical tools that helped define morphogenesis as an emergent property of 
physical and biochemical interactions, as well as cyclical, multidimensional, and 
nonlinear feedback schemes operating at the level of molecules, genes, and proteins, 
and extending up through tissues. In stark contrast to the neo-Darwinian view of 
the relationship between genotype and phenotype, Alberch argued that “genes are 
just one step in the chain of interactions; gene expression is both the cause and the 
effect of a morphogenetic process” (Alberch 1991:6). Using amphibian limb buds 
as a model system for studying the relationship between construction rules and 
morphological outcomes, Alberch and his colleagues experimentally manipulated 
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parameters such as cell number (using the mitotic inhibitor colchicine, for example) 
and showed that changes in size and shape of the limb, and number of the digits 
became altered in a non-random way once a critical threshold was reached (Alberch 
and Gale 1983, 1985; Shubin and Alberch 1986). 

While these studies predated the technical ability to link such outcomes with 
molecular biology (specifically, underlying changes in gene expression), their results 
were completely consistent with predictions made in their mathematical models and 
pattern-generating algorithms (Oster et al. 1988), and showed that the phenotypes 
arising from perturbations to developmental programs were not stochastic. Because 
of such findings, Alberch argued that “even if the parameters of the system are ran- 
domly perturbed, by either genetic mutation or environmental variance or experi- 
mental manipulation during development, the system will generate a limited and 
discrete subset of phenotypes. Thus the realm of possible forms is a property of the 
internal structure of the system” (Alberch 1989:27). Analyses of genetic mutations 
and experimental manipulations in a range of model organisms by many subsequent 
workers in the field provided critical information on the internal structure of devel- 
opmental systems. In particular, these types of approaches have been especially pro- 
ductive with regard to understanding how parameter changes in construction rules 
on the molecular and cellular levels have likely played a generative role during the 
evolution of size and shape in the vertebrate skull. 


7.3 TIME, SIZE, AND SHAPE IN THE VERTEBRATE SKULL 


For numerous reasons, including its inimitable paleontological record, its measurable 
geometry, its evolutionary adaptability, its functional significance, and its easily visu- 
alized embryogenesis, the vertebrate skull has long been the subject of intensive 
research on size and shape (de Beer 1937; Hanken and Hall 1993). This has occurred 
chiefly with regard to; (a) genes that affect skeletal element identity (Balling et al. 
1989; Lufkin et al. 1992; Gendron-Maguire et al. 1993; Rijli et al. 1993; Schilling 
1997; Qiu et al. 1997; Hunt et al. 1998; Smith and Schneider 1998; Pasqualetti et al. 2000; 
Grammatopoulos et al. 2000; Depew et al. 2002; Creuzet et al. 2002; Kimmel et al. 
2005); (b) tissue interactions required for mesenchymal differentiation into car- 
tilage and bone (Schowing 1968; Tyler 1978; Bee and Thorogood 1980; Hall 1980, 
1982b; Tyler 1983; Thorogood et al. 1986; Thorogood 1987; Hall 1987; Richman and 
Tickle 1989, 1992; Dunlop and Hall 1995; Shigetani et al. 2000; Ferguson et al. 2000; 
Couly et al. 2002; Francis-West et al. 2003; Merrill et al. 2008); (c) secreted mol- 
ecules that regulate axial polarity and skeletal outgrowth (Barlow and Francis-West 
1997; Francis-West et al. 1998; Schneider et al. 2001; Hu et al. 2003; Abzhanov and 
Tabin 2004; Crump et al. 2004; Wilson and Tucker 2004; Wu et al. 2004; Abzhanov 
et al. 2004; Liu et al. 2005; Marcucio et al. 2005; Wu et al. 2006); and (d) mesenchymal 
control of species-specific pattern (Andres 1949; Wagner 1959; Noden 1983; Schneider 
and Helms 2003; Tucker and Lumsden 2004; Mitsiadis et al. 2006). 

The special ability of mesenchyme to transmit species-specific information on 
size and shape has been recognized primarily through interspecific grafting experi- 
ments of mesenchymal cells destined to form the jaw skeleton (Noden and Schneider 
2006; Lwigale and Schneider 2008). The exact molecular mechanisms through 
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which mesenchyme performs this complicated function appear to involve the ability 
of mesenchyme to determine the timing of its own gene expression and differentia- 
tion, as well as that of adjacent tissues such as epithelia (Schneider and Helms 2003; 
Eames and Schneider 2005; Schneider 2005; Merrill et al. 2008). Taken together, 
results from genetic, molecular, and cellular studies lead to the conclusion that the 
regulation of skeletal size and shape by mesenchyme involves multiple gene regula- 
tory networks, reciprocal signaling interactions with adjacent tissues, and hierarchi- 
cal levels of control. 


A Biro BEAks 


Studies on the beaks of birds have been particularly helpful in identifying factors that 
influence skeletal size and shape (Helms and Schneider 2003; Schneider 2005, 2007; 
Fish and Schneider 2014c; Schneider 2015). For example, differential domains of 
Bmp4 expression in beak progenitor cells underlie variation in beak depth and width 
among birds including Darwin’s finches, cockatiels, chicks, and ducks (Abzhanov 
et al. 2004; Wu et al. 2004, 2006). 

Beak length seems to be managed separately through a calmodulin-dependent 
pathway (Abzhanov et al. 2006; Schneider 2007). Similarly, factors including 
SHH, FGFs, WNTs, and BMPs, which are secreted from adjacent epithelial tissues 
also appear to affect the shape and outgrowth of the beak skeleton (MacDonald 
et al. 2004; Young et al. 2014; Hu and Marcucio 2009; Foppiano et al. 2007; 
Hu et al. 2015a, b; Hu and Marcucio 2012; Brugmann et al. 2007; Brugmann 
et al. 2010; Abzhanov and Tabin 2004; Bhullar et al. 2015; Grant et al. 2006; 
Wu et al. 2006; Ashique et al. 2002a; Richman et al. 1997; Rowe et al. 1992; 
Szabo-Rogers et al. 2008; Mina et al. 2002; Doufexi and Mina 2008; Havens 
et al. 2008; Schneider et al. 1999, 2001). A clearer picture of how these signaling 
pathways are regulated and how changes to their regulation affect skeletal size 
and shape has begun to emerge. 


B = Quatn—Duck CHIMERAS 


In particular, additional details on molecular and cellular mechanisms through 
which the craniofacial skeleton acquires its proper size and shape have come 
from our studies, using a unique avian chimeric transplantation system that exploits 
species-specific differences between Japanese quail and white Pekin duck (Schneider 
and Helms 2003, 2005, 2007; Lwigale and Schneider 2008; Jheon and Schneider 
2009; Ealba and Schneider 2013; Fish and Schneider 2014b). 

As a proxy for studying the orchestration of morphogenesis more generally, we 
have been posing the question of how do skeletal elements in the jaw skeleton of quail 
and duck achieve their distinct size and shape? Quail have short and narrow jaws in 
comparison to those of duck, which are relatively long and broad (Figure 7.1a, b). We 
have focused on the lower jaw (Figure 7.1c), which forms embryonically within the 
paired mandibular primordia. Neural crest mesenchyme (NCM) that migrates from 
the caudal midbrain and rostral hindbrain is the only source of precursor cells that 
give rise to cartilage and bone within the skeleton of the face and jaws (Figure 7.1d) 
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FIGURE 7.1 The quail—duck chimeric system for studying time, size, and shape in the head 
skeleton. (a) Head skeletons of adult Japanese quail (Coturnix coturnix japonica) and (b) 
white Pekin duck (Anas platyrhyncos) showing species-specific differences in size and shape. 
(c) Lower jaws of adult duck and quail. (d) Neural crest cells generate the facial and jaw 
skeletons (blue) whereas mesoderm forms the caudal cranial vault and skull base (orange). 
(e) Schematic of an embryonic rostral neural tube (dorsal view) depicting the origin of neural 
crest mesenchyme (NCM) from the forebrain (fb), midbrain (mb), and hindbrain (hb). NCM 
destined for the jaw primordia are grafted (green arrow) from a quail donor (red) to a duck host 
(blue). (f) Embryonic quail (red squares) and duck (blue circles) have distinct rates of matura- 
tion but can be stage-matched for surgery (green triangle on Y-axis) by setting eggs in the incu- 
bator at separate times. Approximately three embryonic stages distinguish faster-developing 
quail from duck embryos within two days following surgery, and this three-stage difference 
remains relatively constant during the period of jaw morphogenesis. (Panels [a, b] modified 
from Fish, J.L. et al., Development, 141, 674-684, 2014; Panels [c, e, f] modified from Eames, 
B.F., and Schneider, R.A., Development, 135, 3947-3958, 2008; Panel [d] modified from 
Schneider, R.A., J. Anat., 207, 563-573, 2005, and based on a drawing from D. Noden.) 
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(Le Lièvre and Le Douarin 1975; Noden 1978; Couly et al. 1993; Köntges and 
Lumsden 1996; Helms and Schneider 2003; Noden and Schneider 2006). 

Our experimental strategy involves transplanting pre-migratory NCM between 
quail and duck embryos (Figure 7.le). We transplant NCM either bilaterally, so that 
donor cells fill both sides of the host jaw skeleton, or unilaterally, which allows 
the nonsurgical side of the host to serve as an internal control. Unilateral trans- 
plants enable us to compare donor- and host-derived tissues in the same chimeric 
embryo (Tucker and Lumsden 2004; Eames and Schneider 2005, 2008; Lwigale and 
Schneider 2008; Fish and Schneider 2014b; Solem et al. 2011; Tokita and Schneider 
2009). A powerful and serendipitous feature of this chimeric system is the fact that 
quail embryos develop at a quicker rate than duck embryos (17 vs. 28 days from 
fertilization to hatching), which causes faster-developing quail cells and relatively 
slower-maturing duck cells to interact with one another over time while they become 
progressively asynchronous (Figure 7.1f). Having such divergent developmental tra- 
jectories conveniently offers a way to screen for the effects of donor cells on the 
host by looking for species-specific changes to the timing of gene expression, cell 
differentiation, and tissue formation. Consequently, and especially in the context 
of the aforementioned holy trinity of evolutionary developmental morphology, this 
system also affords us with the unique opportunity to evaluate directly and in the 
same embryo, the effects of changes in growth rates and the timing of develop- 
mental events on size and shape. We can use an anti-quail antibody, which does 
not recognize duck cells, to distinguish the contributions of donor versus host and 
we can quantify the proportion of quail versus duck cells at the molecular level, 
using a PCR-based strategy (Schneider 1999; Lwigale and Schneider 2008; Ealba 
and Schneider 2013; Fish and Schneider 2014b; Fish et al. 2014; Hall et al. 2014; 
Ealba et al. 2015). 

Once quail and duck cells become mixed within chimeras, they become chal- 
lenged to assimilate two separate morphogenetic programs controlling species- 
specific size and shape. Chimeric “quck” contain quail donor NCM inside of a duck 
host whereas “duail” have duck NCM in a quail host. As a result, we can discover 
mechanisms directing jaw patterning by (1) characterizing donor-mediated changes 
to jaw size and shape; (2) assaying for temporal and spatial shifts in developmen- 
tal events underlying cartilage and bone formation such as mesenchymal conden- 
sation and differentiation; (3) analyzing the effects of donor NCM on non-NCM 
host derivatives that participate in skeletal patterning and growth such as epithelia, 
muscles, blood vessels, and osteoclasts; (4) looking for genes that become differen- 
tially expressed in chimeras; and (5) modulating the expression of these genes to test 
if they account for the chimeric phenotype and affect skeletal size and shape (Eames 
and Schneider 2005, 2008; Noden and Schneider 2006; Merrill et al. 2008; Tokita 
and Schneider 2009; Solem et al. 2011; Hall et al. 2014; Ealba et al. 2015). 

An important point to emphasize is that this chimeric system can reveal in a more 
or less “normal” developmental context those molecular and cellular interactions 
between the donor and host that are divergent and ultimately generative of species- 
specific size and shape. In this context, the quail-duck chimeric system offers a 
unique opportunity to observe what Alberch et al. (1979) called “the morphoge- 
netic unfolding of genetic programs in ontogeny and their alteration in the course of 
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phyletic evolution” (p. 297). We can also probe for construction rules and identify 
those parameter changes that may account for evolutionary differences between each 
species. Along similar lines, Shubin and Alberch (1986) argued that while the basic 
morphogenetic rules have remained the same; what have changed during vertebrate 
evolution are the parameters through which these interactions occur (Etxeberria and 
De la Rosa 2009). 

By combining a classical comparative method (Sanford et al. 2002) with experi- 
mental embryology (i.e., the quail—duck chimeric transplant system), we have found 
that NCM relies upon multiple mechanisms to exert cellular control over time, size, 
and shape, primarily through three phases of development: 


e At the onset of NCM migration, quail and duck embryos allocate different 
numbers of progenitors to the presumptive jaw region, with duck having 
significantly more cells (Fish et al. 2014). 

e Thereafter, when these populations of NCM expand, there is species-specific 
regulation of, and response to, critical signaling pathways in a manner that 
is dependent on their own rates of maturation (Eames and Schneider 2008; 
Merrill et al. 2008; Hall et al. 2014). 

e Lastly, as these progenitors start to form cartilage and bone, they execute 
autonomous molecular and cellular programs for matrix deposition and 
resorption through patterns and processes that are inherent to each spe- 
cies and deeply rooted in the timing of developmental events (Eames and 
Schneider 2008; Merrill et al. 2008; Mitgutsch et al. 2011; Hall et al. 2014; 
Ealba et al. 2015; Schneider 2015). 


Moreover, on the molecular level, the SHH, FGF, BMP, and TGF signaling pathways 
all seem to be clearly but not unexpectedly involved since many members and targets 
show species-specific expression and they become altered in quail—duck chimeras. 

Thus, the ability of NCM to regulate the timing, levels, and spatial patterns of 
gene expression and to do so in a species-specific manner, likely modulates the pro- 
liferation, differentiation, and growth of skeletal progenitors, and determines the 
size and shape of cartilage and bone. Such work offers insight into the many ways 
NCM exerts its regulatory abilities during ontogeny and phylogeny, which has been a 
long-standing question in the field (Gans and Northcutt 1983; Noden 1983; Maderson 
1987; Hall and Horstadius 1988; Hanken 1989; Baker and Bronner-Fraser 1997; Hall 
1999, 2000b; Graham 2003; Santagati and Rijli 2003; Trainor et al. 2003; Graham 
et al. 2004; Le Douarin et al. 2004; Noden and Schneider 2006; Jheon and Schneider 
2009; Fish and Schneider 2014c; Sanchez-Villagra et al. 2016). Specific examples of 
the multiple ways NCM exercises control over skeletal size and shape, especially by 
keeping track of time, are detailed in the sections below. 


7.4 EARLY CELLULAR DETERMINANTS OF JAW SIZE AND SHAPE 


The genesis of NCM involves several sequential embryonic events, including induc- 
tion at the boundary between neural and non-neural ectoderm, specification and 
regionalization along the dorsal margins of the neural folds, regulation of cell cycle 
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and maintenance of multipotency, transition from epithelium to mesenchyme, and 
migration throughout the head and trunk (Betancur et al. 2010; Nikitina et al. 2008). 
NCM that arises from the midbrain through the first and second rhombomeres of 
the hindbrain migrates into the mandibular primordia (Le Lièvre and Le Douarin 
1975; Noden 1978; Couly et al. 1993; Kontges and Lumsden 1996; Schneider et al. 
2001). While the gene regulatory networks and developmental programs that govern 
these morphogenetic events are extremely conserved across vertebrates (Nikitina 
et al. 2008; Depew and Olsson 2008; Bronner-Fraser 2008; Northcutt 2005), much 
remains to be understood about when and where changes can lead to the evolution of 
species-specific morphology. 


A — QUANTIFYING JAW PRECURSOR CELLS 


For this reason, we have been concentrating on exactly when and where ducks assem- 
ble their long bills, compared to quail, who make relatively short beaks. Using a 
simple analogy that constructing a taller building might involve adding more bricks, 
as opposed to bigger bricks (Fish and Schneider 2014a), and following the spirit of 
Alberch and his construction rules, we set out to determine the number of jaw pre- 
cursors that migrate into the mandibular primordia in duck versus quail (Figure 7.2a, b). 
We started by counting NCM at key embryonic stages (Fish et al. 2014). 

At an initial embryonic stage, when NCM is specified along the neural folds, duck 
and quail appear to have the same total amount of cranial NCM. But soon afterwards, 
when NCM coalesces on the dorsal margins of the neural tube, duck has about 15% 
more NCM at the midbrain and rostral hindbrain levels, which is where the popula- 
tion that migrates into the presumptive jaw region originates. Remarkably only sev- 
eral stages later, the mandibular primordia of duck contain twice as many cells as do 
those of quail (Figure 7.2c). To understand how a 15% difference could quickly lead 
to a doubling in size, we looked for species-specific variation in cell proliferation and 
cell cycle length. Our results show that although duck has a longer cell cycle (13.5 ver- 
sus 11 hours in quail), if the total duration of each embryonic stage during this period 
is taken into account in terms of absolute time (45 versus 32 hours), then duck cells, 
in fact, proliferate more than those of quail. By maintaining their intrinsic rates of 
maturation, duck deploy a cellular mechanism that increases jaw size progressively 
throughout development (Fish and Schneider 2014c; Fish et al. 2014; Schneider 2015). 
In so doing, they directly link developmental time with size. 


B Genes AND BRAIN REGIONALIZATION 


But how might duck initially generate more midbrain NCM that can then migrate 
into the presumptive jaw region? To address this question, we assayed for species- 
specific variation in the expression of genes known to affect brain regionalization. 
We examined the expression of Pax6 in the forebrain, Otx2 in the forebrain and 
midbrain, Fgf8 at the midbrain—hindbrain boundary, and Krox20 in rhombomeres 3 
and 5 of the hindbrain (Figure 7.2d). Landmark-based morphometrics was used to 
compare brain shape in duck and quail embryos after neurulation, and we pin- 
pointed species-specific differences that were correlated with changes in domains 
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FIGURE 7.2 Molecular and cellular control of species-specific size and shape. (a, b) Frontal 
views of the heads of quail and duck embryos showing differences in the size of the mandibular 
primordia (ma), from which the lower jaw skeleton develops. (c) At this stage, the mandibular 
primordia (ma) in quail embryos (red) is approximately half the size of that of duck (blue) in 
terms of total number of cells. (d) Quail and duck embryos have distinct shapes and region- 
alization of the rostral neural tube. The duck midbrain (mesencephalon) is foreshortened and 
broader mediolaterally. Genes including Foxg/, Pax6, Otx2, Fgf8, and Krox20 are expressed 
in specific domains, each domain being shifted more anteriorly in duck than in quail. The 
lower jaw skeletons of duck (e) and (f) quail show stage-specific and species-specific differ- 
ences in size and shape with duck being longer and more curved. (g) In quck mandibles, the 
quail donor-derived jaw skeleton (red) is shorter and straighter than the contralateral duck 
host-derived jaw skeleton (blue), which is longer and curved. (Panels [a, b, e, f, g] modi- 
fied from Fish, J.L. et al., Development, 141, 674—684, 2014; Panel c modified from Merrill, 
A.E. et al., Development, 135, 1223-1234, 2008; Panel [d] modified from Schneider, R.A., 
Curr. Top Dev. Biol., 115, 271-298, 2015.) 
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of gene expression. Most strikingly, we found that the duck midbrain is shorter and 
broader and has a correspondingly restricted domain of Otx2 expression along the 
anterior to posterior axis. Presumably, a broader midbrain in duck allows more 
NCM to accumulate and migrate into the mandibular primordia. Importantly, 
we observed these differences in Otx2 expression even before neural tube for- 
mation or the genesis of NCM, demonstrating that species-specific patterning 
mechanisms affecting jaw size may function at the earliest stages of development. 
Thus, this work reveals how small spatial and temporal modifications to aspects 
of developmental programs controlling the allocation and proliferation of NCM 
have likely influenced the course of jaw size evolution. 


C REGULATION OF JAW LENGTH 


In addition to discovering that the total amount of NCM present in the mandibular 
primordia is a determinant of species-specific jaw size, we also paradoxically found 
that reducing or augmenting NCM by up to 25% does not significantly alter jaw 
length (Fish et al. 2014). This is consistent with other experiments showing that the 
jaw can return to its normal size after extirpation of precursor cells at the level of 
the neural folds, a process often referred to as regulation (Scherson et al. 1993; Hunt 
et al. 1995; Sechrist et al. 1995; Couly et al. 1996). In these previous investigations, 
however, normal jaw length was thought to arise from regeneration of NCM along 
the neural tube, either by a re-specification of remaining dorsal neuroepithelium 
(Sechrist et al. 1995; Hunt et al. 1995), or by an expansion of NCM generated by 
adjoining neural folds (Scherson et al. 1993; Couly et al. 1996). Instead, we con- 
clude that NCM does not regenerate at the level of the neural tube and therefore, the 
restoration of normal jaw length depends upon another compensatory mechanism 
possibly involving signaling interactions with surrounding epithelia. That is to say, 
normal jaw length may also be affected by local regulation of proliferation within 
the postmigratory environment of the mandibular primordia. 

This type of regulative development in the local environment would allow for 
compensation of deficiencies in NCM up to some intrinsic species-specific popula- 
tion size. Such findings are consistent with prior tissue regeneration and transplanta- 
tion experiments revealing that individual organs possess autonomous determinants 
of size and can regulate growth appropriately in various contexts (Stern and Emlen 
1999; Leevers and McNeill 2005). Moreover, that a strong correlation exists between 
innate rates of growth and overall size is well established in birds (Starck 1989; 
Ricklefs and Starck 1998; Starck and Ricklefs 1998). 


7.5 CELLULAR CONTROL OF JAW SIZE AND SHAPE 
DURING SKELETAL DIFFERENTIATION 


Once appropriate amounts of NCM are allocated to the mandibular primordia of 
quail versus duck, the next question is how these differences are integrated into the 
programs for skeletal differentiation that eventually produces species-specific size 
and shape? To answer this question, we examined the formation of Meckel’s carti- 
lage in the lower jaw skeleton (Eames and Schneider 2008). 
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A Mecket’s CARTILAGE AND SPECIES-SPECIFIC SIZE AND SHAPE 


Meckel’s cartilage is more-or-less a cylindrical rod that is derived exclusively from 
NCM (Helms and Schneider 2003; Noden and Schneider 2006; Noden 1978, 1982; 
Noden and Trainor 2005) and rarely goes on to ossify (Kavumpurath and Hall 1990; 
Ekanayake and Hall 1994; Eames et al. 2004; de Beer 1937). To identify molecu- 
lar and cellular mechanisms through which Meckel’s cartilage acquires its species- 
specific size and shape, we unilaterally transplanted NCM from quail embryos into 
a stage-matched duck. These quail donor NCM filled the right half of the duck host 
mandible, which allowed for an unambiguous comparison of donor quail-derived 
versus host duck-derived Meckel’s cartilage development in the same chimeric 
mandible. 

During normal growth of Meckel’s cartilage, conspicuous stage-specific and 
species-specific differences in size and shape emerge in quail and duck. At early 
embryonic stages in both quail and duck, Meckel’s cartilage goes from being slightly 
curved to more S-shaped. Shortly afterward, however, Meckel’s cartilage in duck 
remains curved (Figure 7.2e), whereas Meckel’s cartilage in quail becomes straight- 
ened (Figure 7.2f). Meckel’s cartilage grows in each successive stage thereafter, 
but steadily gets larger in duck. In quck chimeras, quail donor NCM maintained 
its faster rate of growth within the relatively slower duck host environment, and 
Meckel’s cartilage on the donor side was always accelerated by approximately three 
stages. Moreover, the size and shape of the donor side was consistently more quail- 
like compared to that observed on the contralateral duck host side (Figure 7.2g). 
Using landmark-based morphometrics and a Procrustes analysis (Chapman 1990; 
Rohlf and Bookstein 1990; Coppinger and Schneider 1995; Marcus 1996; Roth and 
Mercer 2000; Schneider and Helms 2003; Zelditch 2004), we quantified changes in 
Meckel’s cartilage and found that NCM controls both stage-specific and species- 
specific size and shape. 

To clarify the molecular and cellular mechanisms through which NCM accom- 
plishes this complex task, we assayed for changes in the program of cartilage differ- 
entiation that might presage the genesis of size and shape. Cartilage differentiation 
involves the condensation of pre-chondrogenic mesenchyme, followed by overt 
chondrification where an abundant extracellular matrix is secreted by chondro- 
cytes (Eames et al. 2003; Hall 2005). In quck chimeras, NCM on the donor side 
differentiated into chondrocytes and formed cartilage, following the timeframe of 
quail. Donor-dependent acceleration to the timing of cartilage differentiation was 
evident from the beginning of mesenchymal condensation. The transcription factor 
Sox9, which is the earliest known molecular marker of chondrogenic condensations 
(Healy et al. 1996; Zhao et al. 1997; Eames et al. 2003, 2004), and Col2a1, which is 
directly regulated by Sox9 (Bell et al. 1997), were expressed prematurely by quail 
donor NCM relative to the duck host. We also observed that FGF signaling, which 
operates upstream of Sox9 and chondrogenesis (Bobick et al. 2007; Govindarajan 
and Overbeek 2006; Murakami et al. 2000; Petiot et al. 2002; Healy et al. 1999; 
de Crombrugghe et al. 2000; Eames et al. 2004) was similarly regulated by donor 
NCM in temporal and spatial patterns like those observed in quail. For example, 
while the secreted ligands Fgf4 and FgfS were expressed continuously by duck 
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host epithelium prior to and during formation of Meckel’s cartilage, the receptor 
Fgfr2 was prematurely expressed just by quail donor NCM. When we inhibited FGF 
signaling during this discrete temporal window of receptor activation, we blocked 
chondrogenesis. Thus by exerting control over the timing of FGF signaling and 
the expression of downstream targets such as Sox9 and Col2al, NCM likely trans- 
mits information establishing stage-specific and species-specific size and shape to 
Meckel’s cartilage. 


B EPITHELIA AND CARTILAGE PATTERNING 


While these experiments demonstrate that NCM dictates the size and shape of carti- 
lage, other studies have shown that adjacent epithelia also play essential roles during 
cartilage pattern formation. For instance, in the 1980s Peter Thorogood advanced 
a “flypaper model” in which he proposed that interactions between epithelia and 
mesenchyme drive the production of extracellular matrix, which adhesively “traps” 
migrating NCM at their site of differentiation and leads to the induction of cartilage 
(Garrod 1986; Thorogood 1988, 1993). In the head, such epithelia are associated with 
the surface ectoderm and pharyngeal endoderm around the facial primordia, as well 
as the brain and sensory capsules, all of which are known to initiate and maintain 
chondrogenesis at one stage or another (Thorogood et al. 1986; Hall 1980, 1981). 
While some data suggest that epithelia can play an inhibitory role during 
chondrogenesis (Mina et al. 1994), additional studies demonstrate that epithe- 
lia impart axial polarity and regional identity to the underlying NCM-derived 
skeletal tissues. More specifically, epithelia around the developing jaws and face 
(e.g., frontonasal, maxillary, mandibular primordia) seem to provide positional 
cues and maintenance factors necessary for patterned outgrowth of individual 
skeletal elements along the proximodistal, mediolateral, and dorsoventral axes 
(Hu et al. 2003; Foppiano et al. 2007; Hu and Marcucio 2009; Schneider et al. 
1999; Young et al. 2000; Cordero et al. 2002; Helms and Schneider 2003; Young 
et al. 2010; Chong et al. 2012; Hu et al. 2015a). Experiments that rotate epithelium 
in the mid- and upper face, for example, cause mirror image duplications of distal 
upper beak structures along the dorsoventral axis (Hu et al. 2003; Marcucio et al. 
2005). Similarly, transplantation studies and genetic analyses demonstrate that 
endodermal epithelium lining the pharynx transmits region-specific polarity and 
segmental identity to NCM, which is critical for the proper growth and orienta- 
tion of bone and cartilage in the jaw skeleton (Couly et al. 2002; Haworth et al. 
2007; Kikuchi et al. 2001; Kimmel et al. 1998; Veitch et al. 1999; Piotrowski and 
Nusslein-Volhard 2000; Miller et al. 2000; David et al. 2002; Crump et al. 2004). 
When either ectodermal or endodermal epithelia are rotated surgically, the under- 
lying NCM-derived skeleton follows accordingly. Taken together, such studies 
indicate that the primary role for epithelia is to contribute local signals for gen- 
eralized anatomical pattern, which in turn induce and/or maintain programmatic 
responses from underlying NCM (Richman and Tickle 1989; Langille and Hall 
1993; Tucker et al. 1999; Ferguson et al. 2000; Mitsiadis et al. 2003; Santagati and 
Rijli 2003; Le Douarin et al. 2004; Wilson and Tucker 2004; Fish and Schneider 
2014c; Foppiano et al. 2007; Hu and Marcucio 2009, 2012; Marcucio et al. 2011; 
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Hu et al. 2015b). Importantly, the timing of expression of these epithelial signals 
is under the regulatory control of NCM (Schneider and Helms 2003; Eames and 
Schneider 2005; Merrill et al. 2008). 

The finding that NCM executes autonomous molecular and histological programs 
for cartilage size and shape can be combined with other experimental results about 
the role of epithelia in the following way. If the steps of skeletal patterning involve 
mesenchymal migration, proliferation, condensation, overt chondrocyte differentia- 
tion, and ultimately the morphogenesis of cartilage as a three-dimensional structure, 
then the interactions with pharyngeal endoderm and facial ectoderm, for example, 
would dictate cartilage orientation and regional identity along the oral cavity. Such 
interactions could happen before mesenchymal condensation and promote and align 
the spatial distribution of pre-chondrogenic mesenchyme. In this context, these epi- 
thelia would be acting instructively initially but then assume a more permissive role 
that facilitates the execution of NCM-dependent programs and enables chondrogen- 
esis to proceed in a time-independent manner. 

So, while epithelia derived from ectoderm and endoderm may define where 
chondrogenic condensations occur along an axis, which is presumably quite similar 
between quail and duck, our transplants reveal that NCM consequently responds via 
intrinsic, stage-specific and species-specific programs that determine cartilage size 
and shape. Equivalent roles have also been postulated for epithelia during osteo- 
genesis of the mandible and other bones (Tyler and Hall 1977; Hall 1978, 1987; 
Bradamante and Hall 1980; Hall 1980, 1981, 1982a, b; Hall and Van Exan 1982; Hall 
et al. 1983; Van Exan and Hall 1984; Merrill et al. 2008). Further support is lent by 
the finding that several chondrogenic signaling pathways including FGFs and BMPs 
are expressed continuously by epithelia prior to and during the arrival of NCM in 
the mandible (Francis-West et al. 1994; Wall and Hogan 1995; Shigetani et al. 2000; 
Mina et al. 2002; Ashique et al. 2002b; Havens et al. 2006; Eames and Schneider 
2008; Merrill et al. 2008). By controlling the timing of receptor activation, in this 
case, for Fgfr2, NCM allows the signal transduction required for chondrogenesis to 
proceed, and by doing so, initiates the program for cartilage size and shape. 

Overall, this integrated perspective on the roles of mesenchyme and epithelium 
in the establishment of size and shape is also consistent with classic embryological 
work from the lab of Hans Spemann who first discovered the origins of species-specific 
pattern in the 1920s and 1930s through interspecific grafting experiments and espe- 
cially by exchanging mouth-forming tissues between frogs and newts (Spemann and 
Mangold 1924; Spemann and Schotté 1932; Spemann 1938; Fassler 1996; Noden 
and Schneider 2006). These remarkable experiments showed that general ana- 
tomical features of the mouth are guided by local signals, but that species-specific 
pattern is dictated by information in the responding cells. Evidently, Spemann 
interpreted his finding to mean that, “The ectoderm says to the inducer, ‘you tell 
me to make a mouth; all right, Pll do so, but I can’t make your kind of mouth; 
I can make my own, and Pll do that”’ (Harrison 1933). Ensuing transplant experi- 
ments between salamanders and frogs, between mice and chicks (in this instance for 
jaws and teeth), as well as divergent species of birds, including quail, chick, duck, 
and emu have also supported the conclusion that species-specific pattern is largely 
driven by NCM (Andres 1949; Wagner 1959; Lumsden 1988; Mitsiadis et al. 2003; 
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Lwigale and Schneider 2008; Sohal 1976; Yamashita and Sohal 1986; Schneider and 
Helms 2003; Tucker and Lumsden 2004; Eames and Schneider 2005; Schneider 
2005; Noden and Schneider 2006; Eames and Schneider 2008; Jheon and Schneider 
2009; Tokita and Schneider 2009; Fish and Schneider 2014c; Fish et al. 2014; Hall 
et al. 2014; Ealba et al. 2015; Schneider 2015). 


C_ Time As A DEVELOPMENTAL MODULE FOR CHONDROGENESIS 


Our results suggest that the program for chondrogenesis through which NCM imple- 
ments species-specific size and shape is integrated at multiple levels and through 
time as a developmental module. This is equivalent to the identification of develop- 
mental modules and the role proposed for mesenchyme in other embryonic systems 
such as epidermal appendages (Eames and Schneider 2005; Schneider 2005). 

In a similar vein as that described by Alberch (Alberch 1982a), modularity is 
predicated on the observation that many developmental programs appear to func- 
tion as semi-autonomous, self-directing, and hierarchical units that can be con- 
tinuously iterated during development and rapidly diversified during evolution as a 
consequence of the inductive relationships among their constituent parts (Raff 1996; 
Bolker 2000; West-Eberhard 2003; Schlosser and Wagner 2004). The notion that 
NCM engineers size and shape by presiding over a highly integrated developmen- 
tal module is supported by the fact that NCM executes autonomous molecular and 
cellular programs for both the formation of cartilage as a tissue and as a three- 
dimensional organ. Importantly, these programs include many of the same gene 
regulatory networks and signaling molecules that operate during NCM specifica- 
tion, proliferation, and differentiation such as members and targets of the BMP and 
FGF pathways, and that affect the size and shape of cartilage in the avian jaw and 
facial skeletons (Francis-West et al. 1994; Wall and Hogan 1995; Mina et al. 1995; 
Ekanayake and Hall 1997; Barlow and Francis-West 1997; Richman et al. 1997; 
Wang et al. 1999; Tucker et al. 1999; Barlow et al. 1999; Shigetani et al. 2000; Mina 
et al. 2002; Ashique et al. 2002b; Abzhanov et al. 2004; Wilson and Tucker 2004; 
Havens et al. 2006; Schneider 2007; Abzhanov and Tabin 2004; Wu et al. 2004, 
2006; Foppiano et al. 2007). 

Further evidence for modularity as a principal mechanism through which NCM 
exerts control over skeletal size and shape relates to the way NCM accounts for time. 
In fact, this is one of the most striking revelations to emerge from the quail—duck 
chimeric system: NCM keeps track of stage-specific and species-specific size and 
shape concurrently. In other words, as quail donor NCM shifts the timing of carti- 
lage differentiation and morphogenesis in the duck to something like that found in 
the quail, these cells generate stage-appropriate and species-appropriate size and 
shape simultaneously. This result offers a novel mechanism that connects skeletal 
development with skeletal evolution vis-a-vis a single population of cells, the cranial 
NCM. Moreover, this melding of time links ontogeny with phylogeny in a manner 
completely consistent with previous theories of heterochrony as a means to under- 
stand transformations in size and shape. 

While historically, heterochrony has been used to describe changes in the tim- 
ing of developmental events between ancestors and descendants (Russell 1916; 
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de Beer 1930), the concept can also concern comparisons of closely related taxa 
(such as quail and duck) and be employed to assess the effects of changes in rates of 
growth on size and shape (Gould 1977; Alberch et al. 1979; Hall 1984; Roth 1984; 
McKinney 1988a; Foster and Kaesler 1988; Klingenberg and Spence 1993; Raff 
1996). This type of growth heterochrony is probably one of many variables intro- 
duced by the faster-developing quail donor NCM in the relatively slower-growing 
duck host. While such an effect would largely arise from intrinsic species-specific 
differences in maturation rates, another variable could be any experimentally 
induced shifts in relative onsets, cessations, and/or durations of molecular and cel- 
lular events during chondrogenesis. Under normal circumstances, these types of 
changes would be considered instances of sequence heterochrony, which is another 
way changes in time can relate to changes in size and shape (Smith 2001, 2002, 
2003), particularly with regard to reciprocal epithelial-mesenchymal interactions 
underlying skeletal evolution (Smith and Hall 1990). 

The predisposition of quail NCM to follow its endogenous rate and time for car- 
tilage development likely arises from cell-autonomous mechanisms that limit the 
cycling and proliferation of cells to a quail-specific timetable. As a result, chondro- 
genesis advances three embryonic stages ahead of schedule, and Meckel’s cartilage 
attains species-specific size and shape. To be clear, this scenario was not the only 
theoretically possible outcome of our transplant experiments. Quail donor NCM 
could have acted naively, followed the timetable of the duck host, and made car- 
tilage that was duck-like in morphology; or they could have become confused and 
created some novel anatomy that was either a combination of, or unlike what is 
normally observed in quail or duck. But instead, within a duck host environment 
and all that entails in terms of duck-specific signaling, quail donor NCM altered the 
relative timing and rates of differentiation, executed an innate program of cartilage 
morphogenesis that replaced and/or superseded the duck program, and in so doing, 
made something like that normally observed in quail. Thus, not only does NCM 
coordinate the developmental timing of its own derivatives, but host epithelium also 
responds to this premature induction and expresses secreted molecules on the donor 
timetable as well (Schneider and Helms 2003; Eames and Schneider 2005; Merrill 
et al. 2008). 

Overall, our work supports the conclusion that heterochrony can underlie the 
species-specific evolution of size and shape, but in the case of the quail—duck chime- 
ric system, such heterochrony does so with at least two important caveats. First, since 
quail donor NCM followed their own timetable and acted as they would normally 
do, the heterochrony we created is not heterochrony in the true sense of the word. 
In this chimeric system, absolute time remained the same, and the heterochrony we 
constructed can only be contemplated in terms of relative timing of developmental 
events (i.e., to that of the duck host). Second, timing is not the only thing that was 
altered in these transplants. Once inside the duck host, quail NCM seemingly and 
progressively implemented a quail-specific genome. Likely, donor NCM does so in 
response to shared common signals present in duck host epithelium (e.g., FGF and 
BMP) that appear to be expressed continuously during a broad developmental win- 
dow, and which might be able to accommodate any difference in stage between the 
donor and the host. Even transplants in avian species with much wider disparities in 
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maturation rates like quail and emu (i.e., 17 vs. 58 days from fertilization to hatch- 
ing), which are separated by approximately seven embryonic stages during chondro- 
genesis, demonstrate that apparently there are few limits in the ability of the host to 
support the deployment of NCM-mediated programs for cartilage and bone at any 
given time during development (Hall et al. 2014). Similarly, duail chimeras, in which 
slower-growing duck NCM act out their programs on a delayed timetable relative 
to faster-developing quail host and consequently make duck-like structures, reveal 
that the same phenomenon is true in reverse (Schneider and Helms 2003; Eames and 
Schneider 2005; Merrill et al. 2008). 


D TIME AS A DEVELOPMENTAL MODULE FOR OSTEOGENESIS 


In a similar manner to what we have observed for cartilage, NCM also appears to 
provide species-specific information on size and shape to the bone in the craniofacial 
skeleton by setting the timing of key events during osteogenesis. Following trans- 
plants of NCM destined to form the lower jaw, quail NCM maintains its faster time- 
table for development and autonomously executes molecular and cellular programs 
that initiate and synchronize each discrete step of osteogenesis including induction, 
proliferation, differentiation, osteoid deposition, mineralization, and matrix remod- 
eling (Merrill et al. 2008; Hall et al. 2014; Ealba et al. 2015). 

Again, this role as a developmental timekeeper holds true both in reverse and in 
the extreme, as evidenced by transplants of duck NCM into quail (i.e., duail) and 
quail NCM into emu (i.e., qumu), respectively. In accordance with one of Alberch’s 
theoretical predictions concerning parameter changes to a construction rule, we find 
that NCM determines the timing of bone formation in the jaw skeleton by controlling 
cell cycle progression. In particular, we observed that NCM regulates the cell cycle 
through stage- and species-specific expression of cyclin and cyclin-dependent kinase 
inhibitors (CKD), including p27 (CdknI1b), which is a CKI that decreases proliferation 
in cell types such as differentiating osteoblasts; cyclin E (Ccne/), which is required 
for G1/S phase transition; and cyclin B1 (Ccnb/), which is required for G2/M phase 
transition (Zavitz and Zipursky 1997; Coats et al. 1996; Drissi et al. 1999). Our 
data suggest that differences in the expression or post-translational processing of 
these cell cycle regulators may enable species such as quail to lessen mesenchymal 
proliferation and form a faster-developing and smaller beak. For example, in quail 
and quck we observed an up-regulation of p27 relative to that observed in duck. 
Previous studies have shown that p27 is correlated with size, including p27-deficient 
mice, which are substantially larger than wild-type littermates and have no appar- 
ent defects in skeletal development (Drissi et al. 1999). Additionally, the developing 
duck frontonasal process has a lower p27 level than in chick (Powder et al. 2012), and 
also the mandibular primordia shows tissue-specific post-translational regulation of 
p27, like what has been reported in other systems (Hirano et al. 2001; Zhang et al. 
2005). Thus, modulating p27 may be a means to influence tissue- and species-specific 
size and/or overall growth. Ultimately, such a direct mechanistic link between the 
regulation of cell cycle progression and the sequence of developmental events during 
osteogenesis likely endows NCM with the capacity to generate changes in skeletal 
size and shape during evolution. 
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Another mechanism through which NCM appears to determine the size and shape 
of bones in the jaw skeleton is through members and targets of TGF and BMP 
pathways, especially, osteogenic transcription factors such as Runx2, which become 
expressed prematurely and at higher levels in quck chimeras (Merrill et al. 2008; 
Ealba and Schneider 2013; Hall et al. 2014; Ealba et al. 2015). Runx2 is a master 
regulator of bone formation since its expression is sufficient to direct osteoblast differ- 
entiation, initiate the timing of mineralization, and affect skeletal size (Eames et al. 
2004; Maeno et al. 2011; Ducy et al. 1997; Komori et al. 1997; Otto et al. 1997; Ducy 
et al. 1999; Pratap et al. 2003; Galindo et al. 2005; Thomas et al. 2004). Moreover, we 
observe premature expression of bone matrix-producing genes such as Collall. 

Consistent with what we observed in our cartilage studies, the systemic environ- 
ment of the duck host seems to be more or less permissive and supports osteogenesis 
independently by supplying circulating minerals and blood vessels. NCM controls 
precisely where and when bone forms by dictating the timing of cell cycle progression 
and by mediating the transition from cell proliferation to cell differentiation. If we 
experimentally induce premature cell cycle exit, we can mimic chimeras by accelerat- 
ing and elevating expression of Runx2 and Colla! (Hall et al. 2014). Experimentally 
increasing and accelerating the timing of Runx2 expression leads to a decrease in 
the size of the beak skeleton like that observed in quail. In effect, this mirrors the 
relationship between endogenous Runx2 levels and species-specific beak size, since 
we also observed higher endogenous expression of Runx2 in quail concomitant with 
their smaller beak skeletons. In fact, by the time the jaw becomes mineralized, Runx2 
levels in quail are more than double those of duck. This supports other studies, which 
have predicted a mechanistic connection between expected Runx2 expression levels 
(based on ratios of tandem repeats in DNA) and facial length in dogs and other mam- 
mals (Fondon and Garner 2004; Sears et al. 2007; Pointer et al. 2012). 

That the timing and levels of Runx2 directly affect the size of the craniofacial 
skeleton, is a finding fulfilling predictions made around 75 years earlier by Huxley, 
Goldschmidt, and de Beer concerning genes that alter the time and rate of develop- 
ment (Huxley 1932; Goldschmidt 1938, 1940; de Beer 1954). Insight into how Runx2 
might play this role comes from in vitro studies in which Runx2 both responds to and 
modulates cell cycle progression through direct and indirect mechanisms, including 
repressing rRNA synthesis, and up-regulating p27 expression (Young et al. 2007; 
Galindo et al. 2005; Thomas et al. 2004; Pratap et al. 2003). 

These studies suggest that NCM controls jaw size by maintaining precise species- 
specific levels of essential transcription factors such as Runx2 and by regulating 
the timing of skeletal cell differentiation. Duck NCM seemingly proliferates more 
slowly and expands in size for longer periods of time before differentiating, which 
then leads to larger skeletal elements. In contrast, quail embryos suppress prolifera- 
tive signals more quickly, exit the cell cycle sooner, and achieve a smaller overall 
beak size. This scenario invokes possible changes to the balance between mesenchy- 
mal proliferation and differentiation during a critical phase of osteogenesis, which 
is condensation (Ettinger and Doljanski 1992; Hall 1980; Hall and Miyake 1992, 
1995). Such changes would likely affect the size, shape, and location of these con- 
densations, which can generate morphological variation in development and evo- 
lution (Atchley and Hall 1991; Dunlop and Hall 1995; Smith and Hall 1990; Hall 
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and Miyake 2000; Smith and Schneider 1998). So, in terms of absolute time, earlier 
osteogenic condensations can lead to smaller skeletal elements and ultimately affect 
their shape through allometric growth. With regard to the differentiation of both 
cartilage and bone, the astonishing ability of NCM to transmit information on size 
and shape across embryonic stages and between species in parallel, lends strong sup- 
port to the notion that development has played a generative role during the course of 
skeletal evolution (Alberch 1982b; Maderson et al. 1982). 


7.6 CELLULAR CONTROL OF JAW SIZE AND SHAPE 
DURING LATE-STAGE GROWTH 


Whereas most of our studies reveal that NCM imparts species-specific size and 
shape to the jaw skeleton by controlling the molecular and cellular programs that 
underlie the induction and deposition of cartilage and bone, we have also found that 
a previously unrecognized but perhaps equally important mechanism affecting size 
and shape is the ability of NCM to mediate the process of bone resorption (Ealba 
et al. 2015). Bone resorption is typically associated with bone deposition and as a 
metabolic process helps maintain homeostasis in the adult skeleton (Filvaroff and 
Derynck 1998; Buckwalter et al. 1996; Hall 2005; Teitelbaum 2000; Teitelbaum et al. 
1997; Nguyen et al. 2013; O’Brien et al. 2008). The extent to which bone resorp- 
tion affects the embryonic skeleton has not been studied extensively except for some 
hypotheses proposing differential fields of resorption to account for changes in size 
and shape that arise during the development of the human jaw skeleton (Enlow et al. 
1975; Moore 1981; Radlanski et al. 2004; Radlanski and Klarkowski 2001). 


A Bone RESORPTION 


Bone resorption relies on the activities of two cell types, which can be distinguished 
by their distinct embryological lineages and morphology: 


Osteoclasts, which come from the mesodermal hematopoietic lineage 
(Jotereau and Le Douarin 1978; Kahn et al. 2009), have traditionally been 
thought of as the principal population of bone-resorbing cells (Hancox 
1949; Martin and Ng 1994; Teitelbaum et al. 1997; Filvaroff and Derynck 
1998; Teitelbaum 2000; Boyle et al. 2003). In our quail-duck chimeric sys- 
tem, all osteoclasts arise entirely from host mesoderm. 

Osteocytes are the second cell type that resorb bone (Belanger 1969; Qing 
et al. 2012; Tang et al. 2012; Xiong et al. 2014; O’Brien et al. 2008; Xiong 
and O’Brien 2012; Akil et al. 2014; Fowler et al. 2017; Jauregui et al. 2016), 
and in the skeleton of the jaws and face form solely from NCM (Helms and 
Schneider 2003; Noden 1978; Le Lièvre 1978). Hence in quail—duck chime- 
ras, osteocytes are derived exclusively from donor NCM. 


Both osteoclasts and osteocytes secrete an enzyme called tartrate-resistant acid phos- 
phatase (TRAP) when they resorb bone (Minkin 1982; Qing et al. 2012; Tang et al. 
2012). Also, osteoclasts and osteocytes express different molecular markers such 
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as Mmp9, which is found in osteoclasts (Reponen et al. 1994; Engsig et al. 2000), 
and Mmp13, which is detected in osteocytes (Johansson et al. 1997; Behonick et al. 
2007; Sasano et al. 2002). Therefore, following the transplant of NCM into the lower 
jaw of chimeric quck, Mmp9 expression would be coming from duck host-derived 
osteoclasts whereas Mmp13 would be expressed by quail donor-derived osteocytes. 

When we examine the initiation of bone resorption in short-beaked quail versus 
long-billed duck we detect significantly higher levels and different spatial domains of 
TRAP, Mmp9, and Mmp/3 in quail, signifying that quail undergo more bone resorp- 
tion than duck, and indicating that elevated resorption may relate to their shorter 
beaks. Correspondingly, chimeric quck have elevated quail-like levels of TRAP, 
Mmp9, and Mmp13 coincident with their quail-like jaw skeletons. This means that in 
chimeric quck, quail donor NCM executes an autonomous species-specific program 
for bone resorption by way of higher Mmp13 expression and TRAP activity, and also 
through upregulation of Mmp9 expression in duck host osteoclasts. This reveals an 
unexpected NCM-mediated mechanism that potentially contributes to the shorter 
jaws of quail and chimeric quck. In other words, levels of bone resorption in bird 
beaks seem to be inversely proportional to jaw size. 

To test if bone resorption is a determinant of jaw size, we used a biochemi- 
cal approach to activate or inhibit resorption by osteocytes and osteoclasts. We 
administered treatments systemically when bone deposition is just starting and 
resorption has not yet begun. Inhibiting resorption causes the quail lower jaw to 
elongate, whereas activating resorption significantly shortens the jaw (Ealba et al. 
2015). Thus, quail and duck express species-specific developmental programs for 
bone resorption that are distinct in terms of levels and spatial domains, these pro- 
grams are governed by NCM, and bone resorption appears to be a contributing 
mechanism establishing beak length. Such experiments point to a novel function for 
NCM-mediated bone resorption, which is to help control species-specific jaw size, 
and they extend previous studies on Darwin’s finches and other species, which argue 
that an important regulator of beak length is the calcium binding protein, calmodulin 
(Abzhanov et al. 2006; Schneider 2007; Gunter et al. 2014). 

This connection is particularly intriguing because calmodulin has been shown 
to regulate osteocytes and osteoclasts in the local environment (Seales et al. 2006; 
Zayzafoon 2006; Choi et al. 2013a, b). Since calcium signaling is known to affect 
bone resorption (Hwang and Putney 2011; Kajiya 2012; Xia and Ferrier 1996; Xiong 
et al. 2014), NCM-mediated bone resorption may serve as another developmental 
mechanism that drives the evolvability of the avian beak more generally (Kirschner 
and Gerhart 1998), and determines jaw size more specifically (Gunter et al. 2014; 
Parsons and Albertson 2009). Additionally, this work suggests that bone resorption 
may act like a rheostat during jaw size evolution and one that is particularly sensi- 
tive to the availability of dietary calcium in the environment, the endocrine effects 
of calcium-dependent hormones, and gradients of calcium signaling within the jaw 
primordia (Schneider 2007). 

The spatial and temporal regulation by the NCM of expression domains for 
genes including Mmp9 and Mmp/3 are likely to affect shape as well, by establish- 
ing local zones of resorption that in effect sculpt the bone and promote or inhibit 
directional growth. Genetic disruptions to these genes and others that affect bone 


192 Cells in Evolutionary Biology 


resorption are known to alter the morphology of the craniofacial skeleton. For exam- 
ple, mice with mutations in Mmp2 have abnormal snouts (Egeblad et al. 2007), and 
defects in jaw morphology are observed in humans with clinical conditions such as 
Spondyloepimetaphyseal dysplasia (i.e., Mmp/3), Juvenile Paget’s disease (i.e., Opg), 
and following treatments with high doses of bisphosphonates, such as zoledronic 
acid, which inhibit bone resorption (Gorlin et al. 1990; Lezot et al. 2014). 

Therefore, the remarkable ability of NCM to exert spatiotemporal control over not 
only the induction, differentiation, deposition, and mineralization of bone (Eames 
and Schneider 2008; Hall et al. 2014; Merrill et al. 2008; Schneider and Helms 
2003), but also the resorption of bone, seamlessly integrates the molecular and cel- 
lular determinants of jaw size and shape, and confers NCM with its unique capacity 
to generate species-specific variation during development and evolution. 


7.7 CONCLUSIONS 


In the beginning of his collected works, Waddington (1975) expressed his “deeply 
ingrained conviction that the evolution of organisms must really be regarded as the 
evolution of developmental systems” (p. 7). A long line of evolutionary developmen- 
tal biologists would certainly concur with Waddington’s viewpoint particularly those 
researchers who have focused on allometry and/or heterochrony as mechanisms to 
explain species-specific transformations in size and shape. A major factor behind 
such transformations clearly involves modifications to the fundamental parameter of 
time, mainly in terms of total developmental time, differential rates of growth, and/ 
or the timing of developmental events. Thus as in good comedy, timing is everything. 

While in many respects Minot, Thompson, Huxley, de Beer, Goldschmidt, 
Waddington, Gould, and Alberch were way ahead of their own time, over the past 
25 years, technological and conceptual advances in genetics, genomics, and molecu- 
lar biology have revolutionized the study of pattern formation during development. 
Many of the genes that regulate basic developmental phenomena have been identi- 
fied, and the processes they guide have been redefined in mechanistic terms. Notably, 
we have come to recognize that the construction rules of embryonic development and 
the genetic and epigenetic architecture required to enforce those rules are shared 
broadly across disparate taxa. This has led to the spread of a common language for 
evolutionary developmental biologists studying the embryos of seemingly diverse 
organisms, including but certainly not limited to mice, chicks, frogs, fish, flies, and 
worms. As a direct result, the pace of research in the field has greatly accelerated 
because discoveries in one species swiftly lead to progress in understanding the 
development of other species. This progress has transformed developmental biology 
from a descriptive science into one that can now explain the complexities of organ 
and tissue development as consequences of known signal transduction pathways and 
transcriptional programs. 

Whether changes to the temporal and spatial programs for development become 
propagated at the level of genomic organization, at the cis-regulatory level of individ- 
ual genes, at the transcriptional and post-transcriptional level through the epigenetic 
activities of non-coding RNA, at the level of nodes within gene regulatory networks, 
at the level of biochemical interactions among gene products such as enzymes and 
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other proteins, at the level of post-translation modification of proteins, at the level 
of diffusion-reaction gradients and thresholds that establish the inductive abilities 
of cells, at the level of cell properties and cell movements, or at the level of physical 
and signaling interactions between tissues, the downstream effects on morphological 
phenotypes can range from subtle to profound. Internal modifications to develop- 
mental programs at any of these hierarchical levels of organization can generate the 
variation necessary for evolution, but they would also be buffered by the robustness 
and stability of the internal networks and nested interactions that ultimately work 
together to generate individual morphological units and ensure fidelity for structural 
and functional integration. 

This feature of developmental systems has allowed the vertebrate craniofa- 
cial complex to be both highly conserved in its basic anatomical organization and 
extraordinarily diversified in its size and shape. Individual morphological units 
within the craniofacial complex can become modified rapidly over time, yet still, 
maintain connections and keep relationships that are required for meeting structural 
and functional demands. By focusing on the molecular and cellular regulation of 
species-specific pattern in the craniofacial complex we hope our work has helped 
pinpoint precisely where and when changes to developmental programs can affect 
the course of morphological evolution. Our experiments in quail and duck embryos 
reveal that NCM plays a special role in generating species-specific pattern in the 
craniofacial complex, by dominating its own signaling interactions with surround- 
ing tissues and by way of autonomous morphogenetic programs that can span and 
accommodate fluctuations in time. Simply because of these virtues, cranial NCM 
has likely endured as a key effectuator of skeletal size and shape during development 
and evolution. 
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